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Abstract 

We study the expansion of the universe at late times in the case that the cosmological 
constant obeys certain scaling laws motivated by renormalisation group running in quantum 
theories. The renormalisation scale is identified with the Hubble scale and the inverse radii of 
the event and particle horizon, respectively. We find de Sitter solutions, power-law expansion 
and super-exponential expansion in addition to future singularities of the Big Rip and Big 
Crunch type. 


1 Introduction 

The recently observed HI El El El El EJ accelerated expansion of the universe has opened new ways 
to think about the connection between cosmology and quantum physics since an entire classical 
origin for the acceleration seems improbable. Quantum theories, on the other hand, potentially 
yield a lot of sources for dark energy (DE), the energy form that is responsible for the cosmic 
acceleration 00121. One of the most prominent candidates for DE is a positive cosmological 
constant (CC) because it appears already in Einstein’s equations of classical general relativity. 
In this context, it can be treated as a perfect fluid with a constant energy density A > 0 and 
an equation of state (EOS) lu = p/ A = —1 that corresponds to a negative pressure p = —A. 
Furthermore, on the quantum level the CC can be interpreted as vacuum energy, however, its 
absolute value has not been predicted reliably in any theory yet, and naive estimates are far too 
large fTU] , 

In this work we investigate the late-time cosmological evolution in a special class of cosmolo¬ 
gies that are based on general relativity with a time-dependent CC and Newton’s constant. 
Classically both constants are truly time-independent, but quantum effects may induce a scale- 
dependence in the form of renormalisation group equations (RGE) for these constants. For 
instance, quantum field theory (QFT) on curved space-times [11,1^ and quantum gravity m 
mm provide a framework for the renormalisation group running. However, the correspond¬ 
ing RGEs are not unique and the identification of the renormalisation scale p, on which they 
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depend, is not fixed by the theories. For phenomenological reasons, we should request that the 
scale /x does not change too much over cosmological time scales. In contrast to DE scenarios 
with a time-dependent EOS mm, here, the EOS of the CC is still exactly u = — 1. However, 
one obtains an effective time-dependent EOS m due to the non-standard scaling of the matter 
energy density with the cosmological time. 

The paper is organised as follows: in Sec. El we describe three scaling laws for the CC. Two of 
them are based on QFT on curved space-times, whereas the third one appears in the context of 
quantum Einstein gravity (QEG) at late cosmological times. In the latter case we also consider 
the running of Newton’s constant. The identification of the renormalisation scale fj, with some 
cosmologically relevant quantity is done in Sec. El We consider three possible choices for fj,: the 
Hubble scale H, the (inverse) radius R of the cosmological event horizon and finally the (inverse) 
radius T of the particle horizon. In Sec. El we discuss how to solve Einstein’s equations on a 
Robertson-Walker background when the CC and Newton’s constant depend on the cosmological 
time. The main part of this work can be found in Sec. El where we investigate the late-time 
evolution of the universe for all combinations of scaling laws and scale choices. Since it is not 
always possible to obtain exact or approximate solutions, some of the given results are numerical. 
Finally, we present our conclusions in Sec. El 


2 Scaling laws 


According to QFT on curved space-times mm the CC and Newton’s constant are subject 
to renormalisation group running like the running of the fine-structure constant in quantum 
electrodynamics. The RGE for the CC, corresponding to the vacuum energy density A, depends 
crucially on the considered quantum fields and their masses. For this kind of running we neglect 
the scale-dependence of Newton’s constant in this work. 


In the following we discuss two different RGEs, the first one follows directly from the 1-loop 
effective action and is given by /x^j = Aq in the MS-scheme, where quantum fields with masses M 
determine the coefficient Aq ~ q^M 4 . Note that the sign of Aq depends on the spins of the fields, 
see, e.g., Ref. m for the exact expression. Upon integration one obtains in the case of constant 
masses 


a (m , M 

= 1 - Qi In—, 

Mo 


M 4 


Qi ~ =p- 


( 1 ) 


Ao Mo A o 

where Ao denotes the vacuum energy density today, when the renormalisation scale /x has the 
value /xo- The scaling law 0 has a big disadvantage since it has been derived in a renormalisation 
scheme that is usually associated with the high-energy regime EDI- It is therefore not known to 
what extend it can be applied at late times in cosmology. In addition, to avoid conflicts with 
observations, the field content has to be fine-tuned to obtain |gi| < 0(1), which we assume for 
the rest of this work. 


The second scaling law follows from a RGE that shows a decoupling behaviour. Considering 
only the most dominant terms at low energy, the corresponding /3-function for A is given by 
M^ = Ai/x 2 , where A 1 ±M 2 is set by the masses M and the spins of the fields. RGEs of this 
kind have been studied extensively in the literature E2- Assuming constant masses and mo to 
be of the order of today’s Hubble scale Hq one obtains 


a (m) 
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where Lq := 1 — L\ and Mp denotes the Planck mass today. Here, the running of A is suppressed 
since \L\ \ <C 1 for sub-Planckian masses M. 

The last scaling laws come from the RGEs in QEG main. In this framework the effective 
gravitational action becomes dependent on a renormalisation scale //. which leads to RGEs for A 
and G. An interesting feature is the occurrence of an ultraviolet (UV) fixed-point [21 hr the 
renormalisation group flow of the dimensionless quantities 1 A/r~ 4 and G/i 2 at very early times 
in cosmology corresponding to /r —► oo. Motivated by strong infrared (IR) effects in quantum 
gravity one has proposed that there might also exist an infrared fixed-point in QEG 1221 leading 
to significant changes in cosmology at late times 2 , where /i —► 0. If this were true, one would 
obtain the RGEs 

A = A = M /or 

Ao /A Go /r 2 ’ 1 J 

where corresponds to A = Ao and G = Go- In the epoch between the UV and IR fixed- 
points, A and G vary very slowly with fi. and we will treat them as constants in this region. 
Therefore, we assume that the scaling laws m are valid from today until the end of the universe. 


3 Choice of the scale 


In order to study the effects on the cosmological expansion due to the scaling laws of Sec. El 
we have to define the physical meaning of the renormalisation scale /i. Unfortunately, the 
theories underlying the scaling laws do not determine the scale explicitly 3 , apart from the usual 
interpretations as an (IR) cutoff or a scale characterising the physical environment (temperature, 
external momenta). In our cosmological setting we will investigate three different choices for /j, 
given by the Hubble scale H , the inverse radius RT 1 of the cosmological event horizon, and 
finally the inverse radius T _1 of the particle horizon. 

For the purposes of this work we consider a spatially flat Robertson-Walker metric given by the 
line element 

ds 2 = df 2 — a 2 (t)dx 2 , (4) 

where the scale factor a(t) depends only on the cosmological time t and not on the spatial 
coordinates x. While the Hubble scale 


m .= 


a(t) 
a(t ) 


( 5 ) 


describes the actual expansion rate of the universe, the horizon scales R and T describe the 
cosmological evolution of the future and the past, respectively. The radius R of the cosmological 
event horizon corresponds to the proper distance that a (light) signal can travel when it is 
emitted by a comoving observer at the time t: 

f°° d t' 

m ■■= «(t) / —r v (6) 

Jt a{t') 


In the case that the universe comes to an end within finite time the upper limit of the integral has 
to be replaced by this time. Similar to the event horizon of a black hole, the cosmological event 


1 In this work A denotes the vacuum energy density corresponding to a CC, whereas in many articles about 
QEG A is the CC itself. 

2 However, it was argued in Ref. 1231 that the region of strong IR effects can never be reached. 

3 In the framework of Refs. m, Newton’s constant was found to depend explicitly on the Hubble scale. 
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Figure 1: In a spatially flat ACDM universe with constant A and G the radii of the event horizon R and 
the Hubble horizon H^ 1 approach the same constant value for t —> oo. At early times the 
particle horizon radius T grows linearly with time and exponentially with time for f —> oo. At 
the time t 0 = 0, 9977,7 1 (today) the relative vacuum energy is given by Rao = 0, 73. 


horizon exhibits thermodynamical properties like the emission of radiation with the Gibbons- 
Hawking temperature m- The analogue to R is the particle horizon radius T, which is given 
by the proper distance that a signal has travelled since the beginning of the world (t = 0): 

G (lb 

T(t):=a(t) -- (7) 

Jo ait') 

In a simple cosmological model, where the universe begins at t = 0 and then evolves according 
to the scale factor a(t ) oc t n with 0 < n < 1, the particle horizon radius reads T(t) = 7/(1 — n ). 
Now we assume that at some time to a de Sitter phase sets in, which corresponds to a scale 
factor o(t) oc ex.p(Ho(t— to)) for t > to- This leads to a radius function T that grows exponentially 
with f at late times: 

T(t > t 0 ) = ^T(to) + exp (H 0 (t - t 0 )) - (8) 

Since the asymptotic behaviour of the scales H, R and T plays a major role in this work, it is 
plotted in Fig. |T|for a universe with dust-like matter and A and G being positive and constant. 
Note that if one of the horizon radii diverges, one says that the horizon does not exist. The 
geometric meaning of the scales given here has been discussed, e.g., in Ref. m , and some 
applications in the context of DE can be found in Refs. m- 


4 Time-dependent A and G 

Applying the scaling laws of Sec. El together with the identification of the scale /x with one of 
the cosmological scales of Sec. 01 we obtain time-dependent functions A(f) and G(t). Einstein’s 
equations G a @ = 8irG(Ag a P + T°^ ) and the Einstein tensor G a & remain unmodified apart from 
the additional time-dependence in A and G. For the late-time epoch we consider the spatially flat 
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metric g a ^ following from Eq. (l4|l . and the energy-momentum tensor T af3 describes a perfect fluid 
with energy density p and pressure p. Additionally, we introduce the variable Q := (1 + 3cu)/2 
which is fixed by the EOS co = p/p of the fluid. Finally, the evolution of the scale factor a[t) is 
given by Friedmann’s equations 


a 2 87t , a 8vr a 

2 ~ + p), - — — G(A — Qp ), 

a z 6 a 6 


which have to agree with the Bianchi identities 

0 = G a0 . 0 = [GAg a P + GT a %, 


( 9 ) 


that imply 

G(A + p) + G(A + p + 2 H(Q + 1 )p) = 0. 

For constant A and G the last equation can be integrated to yield the usual scaling law for the 
matter energy density p oc a -2 ^ +1 ). In general, this is not possible when A and G depend on 
time due to the energy transfer between A and p. It therefore implies an effective interaction 
between the gravitational sector and matter, which is not part of the original Lagrangian. In 
this sense it should be compared with gravitational particle production EH! resulting also from 
the interplay of gravity with quantum physics. 

To determine the scale factor a(t) we combine both Friedmann’s equations © resulting in 


K 0 F(t) 


mm 

Go Aq 


with F(t) := —b Q 
a 


H + (Q + 1)H 2 , 


( 10 ) 


where the constant Kq := [Hao(Q + 1)-Hq] 1 fixed by today’s values of the Hubble scale Ho, 
the relative vacuum energy density Hao and the matter EOS Q. 


5 Late-time evolution 

In this section we study the late-time evolution of the universe with variable A and G , thereby 
assuming from today on the validity of the scaling laws m-m and the correct identification of 
the renormalisation scale p with the scales ©-©■ Our aim is to determine in all nine cases 
the possible final states of the universe. This depends, of course, on the choice of parameters, 
but we restrict ourselves to parameter values which comply vaguely to current observations. 
Today, at the cosmological time to = 0, 99HQ 1 we fix the initial values Ho and Hao = 0, 73 by 
observations. Furthermore, the initial value To of the particle horizon radius would be fixed if the 
past cosmological evolution was known from the Big Bang on. In contrast to this, the value Ro 
of today’s event horizon radius depends on the future cosmological evolution and is treated here 
as a free parameter. In the following we derive some properties of the solutions analytically, in 
particular, we study the stability of (asymptotic) de Sitter solutions and the occurrence of future 
singularities ESI, where the scale factor or one of its derivatives diverge within finite time. For 
simplicity we denote a Big Rip or a Big Crunch by the lowest order divergent derivative that is 
positive or negative, respectively. Since some combinations of scaling laws and renormalisation 
scales lead to complicated equations we derive in some cases only approximate or numerical 
statements. 

In order to solve Eq. m numerically, we have to remove the integrals in the definitions ©, 0 
of R and T by a differentiation with respect to t. Using the relations 

R = RH - 1, T = TH + 1 
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Figure 2: Lambert’s W-function is the inverse function of z = xe x . For — 1/e < z < 0 it has two 
real-valued brunches Wo(z) and W-i(z), and for z > 0 only Wq{z) is real. 


an ordinary differential equation for the scale factor can be obtained and numerically integrated. 
Afterwards we have to check whether the functions R and T, calculated from the numerical 
solution of a(t), agree with R and T that follow directly from Eq. (1101) . If they do not match, 
the numerical solution has to be discarded. Furthermore, solutions involving a negative matter 
energy density are questionable on physical grounds. This happens when the vacuum energy 
density A becomes greater than the critical energy density p c = 3H 2 /(8ttG), which follows from 
the first Friedmann equation ©• 

5.1 A = A 0 (l — <?i In ^), n — H 

For the given scaling law and scale choice Eq. m reads 


K 0 (H + (Q+ 1 )H 2 ) = 1 - qi In 

R o 


( 11 ) 


We first look for asymptotic de Sitter solutions by applying H = 0 and H —> H c . Thus the final 
Hubble scale H 0 is given by 


m 

Hi 


y^AO W u 



( 12 ) 


Here, W u (z) with u = 0,-1 denotes one of the two real-valued branches of Lambert’s W-function, 
which is the solution of z = xe x , see Fig. El For q\ > 0 there is always one solution for H e , 
and for q\ < 0 two solutions exist if the argument of W u in Eq. m is greater than — e 1 . This 
means either q\ < n\ or q\ > ri 2 with n\ := 2/Wo(—H aoC - 1 ) and ri 2 := 2/VF_i(—H aoC -1 ) > n\. 
These solution are stable if 


= -^--2H e K 0 (Q + l) <0, 

where we used Eq. m- For positive q\ this condition is always fulfilled, whereas for negative q\ 
it means llw( fjl o A0 e 2 / qi ) < —1, which follows from Eq. ( 1121 ) . Again, the argument of W u is 


d(K 0 H) 
d H 
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Figure 3: The phase space relation m is shown for q\ = —6; —3, where the universe ends in a Big 
Crunch singularity, and respectively for q\ = — 1;+^, which leads to a de Sitter final state 
( H = 0). In all cases the evolution begins at H = H 0 , where Rao = 0, 73 and H 0 < 0. 


constrained yielding q\ > ri 2 - Therefore only q\ > ri 2 , which includes q\ > 0, leads to stable 
de Sitter solutions. Using the phase space relation ED. which is plotted in Fig. [3 we conclude 
that for other values of q\ the cosmological evolution will always end within finite time in a Big 
Crunch singularity, where H —> 0 and H —> — 00 . 


5.2 A — A 0 (L 0 + ii^), l-i — H 


Mo ' 


Using this choice for A and /r Eq. ED becomes 


— + H\Q — L\(Q + l)fi A0 ) — + 1), 


which has the exact solution 

a(t ) 


do 


sinh([Q + 1][Q A0 L 0 (1 - Ti^ao)] 2 ^) 


with 


n := \{Q + 1)(1 — Lif2 A0 )] -1 . 

At early times, Hqt <C 1, it describes a power-law expansion a(t) oc t n , whereas at late times it 
approaches the de Sitter expansion law a(t) oc exp (H e t) with the final Hubble scale H e given by 


= //, 


Hao-Lo 


1 — LiHao 

Further aspects of this case have been studied, e.g., in Refs. El- 


5.3 AG = A 0 Go^, m = H 

Mo 

Following Sec. El we assume that this scaling law is valid from today on. Here, Eq. ED is given 
by ^ + BH 2 = 0 with B := Q — (Q + 1)Hao- It has an exact power-law solution 

= [(1 + B)(t — fi)] 1 ^' B+1 \ t\ = const., 

«o 

that exhibits a constant acceleration aa/a 2 = — B > 0 as well as = const. A similar solution 
was found in Refs. 123- 
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5.4 A — A 0 (l — <?i In ^), n — R 1 


Since this case has been discussed explicitly in Ref. E2I, we skip some details in the following 
discussion. At first we look for de Sitter solutions, where the inverse Hubble scale and the event 
horizon radius approach R —> R e = H^ 1 in addition to H = 0. Then KqF = 1 + q\ In can be 
solved for R e yielding 


Rl = 


H, 




(13) 


For positive q\ there is always a solution for R e , but it is unstable as can be seen from the following 
consideration. We first insert R = Ro exp([/io-F — l]/gi) and H = \JKqF/[Kq(Q + 1)], which is 
valid in the de Sitter limit, into KqF = qi(H — R _1 ). Then we compute 


d(K 0 F) 

d(K 0 F) 


dS 


^[KqF(Q + 1)] 2 +-i- exp 


1 - K 0 F 

qi 


(14) 


which is positive for q\ > 0 signalling instability. Thus a Big Crunch singularity (F, Ft —> —oo, 
R—y 0 ) might happen for certain values of initial conditions and parameters, or there exists no 
solution. A Big Rip singularity (F, H —> oo, R —> 0) does not exist because it is not compatible 
with R —y oo. 


For negative q\ there are two solutions for R e if the argument of W u in Eq. m is greater 
than —e _1 corresponding to the condition 


Ro > 


2 exp( 2 /qi + 1 ) 

-qiH$Q,Ao 


According to Eq. m it turns out that R e ± := i? e |«=o represents a stable solution and R e 2 : = 
i? e | u= _i an unstable one, respectively. A closer inspection of the solutions R e i /2 shows that R e \ > 
Ro > R e 2 if RqHq£1\o > 1. In the case RqHqQ\o < 1 the parameter range 

2 

Ql < Wq{—\/ (e 1 ^)) 

results to Rq > R e 1 > R e 2 , whereas 


Qi > 


2 

Wo(-l/(e 1 g 3 )) 


leads to R e 1 > R e 2 > Ro, respectively. Moreover, all stable de Sitter solutions are realised except 
the last one corresponding to the parameter region (RqHqQ\o < 1, qi > 2/Wo[— 1 /(e 1 ^ 73 )])- In 
this last case the universe will end in a Big Rip singularity, where F,H—y 00 and R — y 0. A Big 
Crunch (F —> — 00 ) is not possible since it contradicts R —y 00 . 


5.5 A — A 0 (L 0 + /x — R 1 


In the de Sitter limit, H = 0 and R —> R e = H c 1 , we find the solution 

Rl = - L.RD/Lo, 

but it is unstable for all possible values of L\, which follows from 


d(K 0 F ) 

d(KoF) 


dS 


R c HoL 0 n AO > 0 . 
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(a) (b) 

Figure 4: Corresponding to Sec. 15.51 we plot (a) the function K 0 F = L 0 —for L i = ±0, 3 and (b) 
the Hubble scale H{t ) for L\ = +0,1; Rq = 1,1; Hao = 0, 73 and Q = \ (dust), where a Big 
Rip singularity occurs. 


For Li < 0 our numerical calculations did not yield any solution that is compatible with Eqs. (0 
and ED. This is also true in a certain parameter range for L\ > 0, where otherwise a Big Rip 
singularity occurs, where F, H —> oo and R —> 0, see Fig. 0 A Big Crunch is not possible 
since KqF is bounded from below. 


5.6 AG - A 0 G 0 4, M - R- 1 


Here, one cannot find a prediction for R e in the de Sitter limit (H = 0, R —» R e = H^ 1 ) since 
Eq. ED only leads to the constraint Rq = 1/(HqQ,\q). To find solutions corresponding to other 
values of Ro we impose the ansatz a(t) = ao(t — t\) n with ao,n,t\ = const. The Hubble scale 
is given by H = n(t — H) -1 which implies due to Ho > 0 that t\ < t for n > 0 and t\ > t 
for n < 0, respectively. For 0 < n < 1 the event horizon does not exist, in the other cases R can 
be calculated exactly via Eq. ©: 


r f E ,\f 

R{t) = a(t) / — 
Jt a K v 


0 t-hT 


(■ t E - H) (1 ~ n) - (t - fi) (1 ~ n) 
1 — n 


Note that for negative n there is a Big Rip singularity at t t\, where a(t) diverges and the 
universe ends: tE = t\. For n > 1 the scale factor a(t ) describes power-law acceleration and 
there is no future singularity which means tE ~^> oo. As a result, both cases yield 


m 


(t-t i) 
n — 1 


Thus Eq. (HOI) can be solved exactly, which determines the constant n, see Fig. |bj 


ni,2{x) 

with x 


(x-\)± ^J(x - \) 2 - x(x - (Q + 1 )) 
x - (Q + 1 ) 

^ = RlHfa A o(Q + l). 

Aq 


(15) 

(16) 
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Figure 5: The exponents ni^{x) fromEq. 03, which occur in the power-law expansion scale factor a(t) = 
ao(t — t\) n , are shown. In Sec. I5.fi! we find n = n\( Rq/K 0 ), whereas n = ti 2 {Tq/K 0 ) in Sec. 15.91 

For Q > 0 we find (Q + l)” 1 < ri 2 < 1 for all positive values of x, therefore ni can be dropped as 
the event horizon does not exist. In the case Rq < 1/(HqQ^q) the constant n = n\ is negative 
leading to a Big Rip singularity at the time t = t\, and Rq > 1/{HqQ.^q) implies a positive n 
corresponding to power-law acceleration. 


5.7 A = A 0 (l-qiln^), M = T- 1 


Solving Eq. m for this choice of A and fj, can be done approximately. First we propose the 
ansatz a(t) = aoexp(w 2 (t — fi) 2 ) for the scale factor with the constants ao, u and t\. Therefore 


H = 2 u 2 (t — t\), — = 4 u A (t — ti) 2 + 2 u 2 , F = 2 u 2 + 4 u A (Q + l)(t — tp) 2 , 
a 

and the particle horizon radius at late times reads 

T = exp(u 2 (t — ti) 2 ) f exp(— u 2 (t' — ti) 2 )dt r , t > t x = const. 

Jt x 

|^erf (u(t - h)) . 

tx 

For t —> oo the series expansion of the square brackets reads C + 0(t~ l ) with C > 0 and 
thus T ~ C'exp(u 2 (t — ti) 2 ). In this limit the given ansatz for a(t) solves KqF = 1 — giln^ 
exactly and determines the constant u 2 = ^q\F[^Q\o. Therefore we have found an approximate 
late-time solution in the case q\ > 0. 


= exp(u 2 (t — ti) 2 ) • 


The numerical solutions in the case q\ < 0 exhibits a future Big Crunch singularity. At times 
sufficiently before this event the scale factor can be well approximated by the ansatz a(t) = 
aoexp(— u 2 (t — fi) 2 ), which implies F = —2 u 2 + 4u 4 (Q + 1 )(t — t\) 2 and 

t 

, t > t x = const. 

tx 


T = exp(— u 2 (t — ti) 2 ) 


—erf (iu{t - ti)) 
Liu 


We expand the square bracket term around t —> t±, where the scale factor is maximal, and 
find T « C exp(— u 2 {t — fi) 2 ) with C > 0. Solving KqF = 1 — q\ In leads to u 2 = — \qiH^Cl\o 
in accordance with q\ < 0. The numerical and approximate solutions are shown in Fig. |HI 
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Figure 6: Here, the numerical and the approximate solutions for the scale factor a{t) and the Hubble 
scale H(t) from Sec. 15.71 are shown for negative qi = —0, 2. Furthermore, we used Hao = 0, 73; 
TqHq = 3 and Q = \ (dust). In this case a Big Crunch singularity will occur in the future. 
Note that the analytical approximation works very well except near the future singularity. 


5.8 A = A 0 (L 0 + L^), /z = T- 1 


In the late-time de Sitter limit, a{t) oc exp (H e t), the particle horizon radius has the asymptotic 
form 

T oc exp (H e t) —> oo for t —> oo, 

which follows from Eq. This means that KqF —> Lq and thus determines = HqQ\qLq. 
To verify the stability of this solution we calculate 


KqF = -2L, 


T Tl 

J 1 y2 


= -2 L 


n 

2^2 


H + 


At the initial time (today) KqF is negative for L i > 0 and K$F keeps on approaching the 
de Sitter limit unless a sign change in KqF would occur. But this requires a negative F[, which 
is only possible if KqF = H + (Q + 1 )H 2 < 0 at some time. Since KqF > Lq > 0 we conclude 
that the de Sitter limit will be reached always. For L\ < 0 the initial values of KqF and KqF 
are positive and KqF approaches the de Sitter limit. Again, a sign change in KqF is only 
possible for H < o, which requires KqF < 0 at some time. These conditions cannot be realised 
because KqF > 0 as long as KqF > 0. In summary all reasonable values of L\ lead to a stable 
de Sitter final state. 


5.9 AG = A 0 Go^, /* = T 

Mo 


-1 


By using at late times the ansatz a(t) = ao(t — t\) n with a,Q,t,n = const., we find 

F = ( n 2 (Q + 1) - n)(t - fi)~ 2 , 
and the particle horizon radius is given by 


T = a 0 (t - ti) r 


f t0 dt' f l d t' 
Uo a(t') + Jt 0 a(t’)_ 


t —1\ 
1 — n 


+ {t~ ti) r 


Tq- 


(tp ~ h) 
1 — n 


1 —n 


(17) 
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where To = T(to) > 0 depends on the past evolution of the scale factor. In the case n < 0 
a positive Hubble scale H = n(t — ii) _1 requires t, to < t\. which implies T < 0 and thus the 
non-existence of the particle horizon. For 0 < n < 1 the second term in Eq. d can be neglected 
at late times, and Eq. m can be solved exactly, leading to 

0 - \) ± \j(x- ~ x ( x ~(Q + !)) 

x — (Q + 1) 

^L=T$H 2 0 n A0 (Q+ 1). 

To 

This is the same equation as in Sec. l5.6l with Ro replaced by To, see also Fig. El Here, solution ri 2 
is the physical one because of (Q + l)~ x < 77,2 < 1, that corresponds to a decelerating universe 
for t —► 00 . In comparison with Sec. 15.61 the identification of the particle horizon radius as 
renormalisation scale yields a complementary cosmological behaviour. 


ni,2(x) = 

with x := 


6 Conclusions 

In this work we have studied the effect of a variable CC on the late-tinre evolution of the universe. 
Thereby the variability of the CC originates from RGEs that are provided by certain quantum 
theories, leading to a scale-dependence of the CC. By identifying the renormalisation scale with 
cosmologically relevant scales, the CC becomes effectively time-dependent. In contrast to many 
other models of DE, here, the EOS is still that of a classical CC, which implies a non-trivial 
scaling behaviour of the cosmological matter content. Concretely, we have investigated three 
RGEs for the CC and three choices for the renormalisation scale: the Hubble scale and the inverse 
radii of the cosmological event horizon and particle horizon. For one RGE we have also included 
the renormalisation group running of Newton’s constant. In all cases we have determined the 
possible final states of the universe and discussed their properties. Depending on the choice of 
the initial conditions and parameters, we have found at late times asymptotic de Sitter solutions, 
accelerating and decelerating power-law solutions and super-exponential expansion laws of the 
type a(t) ~ exp(cf 2 ). Additionally, we have encountered future singularities of the Big Rip and 
Big Crunch type. However, some solutions should be treated with caution, since they imply 
physically questionable results like negative energy densities. Also, the cosmological behaviour 
near future singularities might get altered by additional effects, that have not been taken into 
account in this work. It was shown, for instance, in Refs. m that quantum effects are able to 
prevent or weaken future singularities. Moreover, higher-derivative terms in the gravitational 
action and a time-dependent Newton’s constant could potentially change the solutions. 

Nevertheless, we stay in this analysis on the basic level and use the results we found as an 
indicator for the (in)stability of the cosmological fate. Indeed, we have found regular solutions 
in many cases, which can be considered as realisable in nature. In this sense, we have found 
de Sitter solutions in Secs. I5~n lA2l IQ and Ol Also the power-law solutions can be called 
regular, they appear for all cases with the scaling law m in Secs. 15.3115.31 and EU Moreover, 
we have found super-exponential expansion laws in Sec. 15.71 Big Rip type future singularities 
exist for all cases with the event horizon as renormalisation scale, see Secs. E21EZD and 15.61 Big 
Crunch solutions, on the other hand, occur only for the scaling law 0 as described in Secs. E3J 
15.41 and 15.71 Finally, a comprehensive overview of our results is given in Tab. Q Since some 
solution types occur only for certain combinations of scaling laws and scales, our analysis helps 
to discriminate between the different cases and discover the nature of DE. In a future work, we 
would like to discuss the consequences of a running CC for the early universe, too. 
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A = A 0 (l - qi ln^h) 

A — Ao(Lo + L\ 

Mo 

A G — 

Mo 

H = H 

dS, BC 

dS 

P 

II 

i 

dS, BR, BC 

BR 

P, BR 

II 

4 

1 

exp(±f 2 ), BC 

dS 

P 


Table 1: This table shows all types of cosmological final states that we have found in Sec. |3 for the all 
combinations of scaling laws for the CC A (and G for the third law) from Sec. |3 and the choices 
for the renormalisation scale fi of Sec. 0 Here, dS denotes de Sitter solutions, P power-law 
solutions, BR Big Rip and BC Big Crunch future singularities, respectively. With exp(±f 2 ) we 
mean the late-tinre behaviour of the scale factor, where a <x exp(— t 2 ) finally leads to a Big 
Crunch. 
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